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Results  of  a study  of  nonequilibrium  effects  in  jet  and  rocket 
nozzles  and  exhaust  plumes  are  presented.  The  study  considered  turbu- 
lence modeling  of  axisymmetric  compressible  jets,  particle  laacn  jets 
and  chemical  nonequilibrium  in  axisymmetric  jets.  It  was  found  tfiat. 
available  turbulence  models  do  not  correlate  a wide  range  of  data  for 
compressible  jets.  Studies  of  particle  laden  jets  provided  reasonabie 
correlacion  to  experimental  data.  The  chemical  nonequilibrium  studies 
suggest  a possible  explanation  for  differences  between  observed  and 
previously  predicted  radiation  behavior  of  axisymmetric  exhaust  plumes. 
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1. 

I . Introduction 

The  fundamental  understanding  of  exhaust  plume  characteristics  is 
of  considerable  practical  importance  to  the  United  States  Air  Forc.e. 

Vehicle  detectability  depends  on  exhaust  plume  radiation,  plume  radar 
cross  sections  affect  vulnerability  of  tactical  weapons,  and  highly 
expanded  afterburning  plumes  can  lead  to  vehicle  contamination  efforts. 

The  detailed  study  of  these  phenomena  with  full  scale  vehicles,  over  the 
range  of  vehicle  operation  is  a difficult  experimental  task.  The  devel- 
opment of  plum.e  prediction  models  provides  an  alternative  means  of  de- 
termining these  phenomena. 

A number  of  codes  have  been  developed  (1,  2,3,4)  which  include 
many  of  the  physical  and  chemical  phenomena  in  the  exhaust  plume.  These 
models  of  exhaust  flow  may  be  divided  into  four  categories  (4);  i)  parallel 
mixing  models,  ii)  uncoupled,  underexpanded  flow  models,  iii)  limited 
coupling  models  and  iv)  complete  coupling  models.  Under  the  assumptions 
of  the  parallel  mixing  model  radial  pressure  gradients  are  considered 
negligibly  small  compared  to  gradients  due  to  turbulent  mixing  and  chem- 
ical reactions.  The  governing  equations  then  take  the  form  of  parabolic 
partial  differential  equations.  In  the  uncoupled  underexpanded  flow 
models  the  flow  is  described  initially  as  inviscid  and  hyperbolic  and  a 
method  of  characteristics  solution  is  employed.  Additional  calculations 
are  necessary  to  assess  viscous  and  chemical  kinetic  effects  which  can  be 
patched  to  the  above  solution.  Limited  coupling  models  allow  the  inclu- 
sion of  coupling  between  some  of  the  governing  differential  equations 
while  neglecting  coupling  between  others.  The  most  general  solution 
represented  by  complete  coupling  has  recently  been  developed  (3).  These 
models  all  rely  on  separate  evaluation  of  chemical  kinetic  and  transport 
property  (oeffirients  as  input  data. 

k.  
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The  study  which  is  reported  here  is  part  of  an  attempt  to  eval- 
uate input  transport  phenomena  and  chemical  effects  in  afterburning 
exhaust  plumes.  The  first  portion  of  this  study  focused  on  the  review 
and  an  attempt  to  evaluate  turbulence  models  of  compressible  axisym- 
metric  homogeneous  jets.  Since  it  is  known  that  particles  are  important 
in  exhaust  pl'jmes,  the  extension  of  one  of  these  models  to  include 
particles  in  the  flow  was  the  next  area  of  the  study.  Chemical  effects 
were  then  added  to  the  mixing  turbulent  jet  flow  to  evaluate  the  effect 
of  initial  mixture  composition  on  the  chemical  history  of  the  jet. 

These  results  are  summarized  briefly  below  before  presenting  the  study 
in  more  detail  in  succeeding  sections  of  the  report. 

Previous  experimental  and  theoretical  studies  of  axisymmetric  jet 
mixing  have  included  compressible  and  incompressible  jets  an  I homogeneous 
and  two  phase  jets.  Homogeneous  incompressible  jets  have  been  stufJied 
most  extensively  and  correlation  between  experiment  and  analyticaj 
descriptions  seems  adequate.  Homogeneous  compressible  jets  have  seen 
considered  (6,  7,8)  and  some  experimental  data  are  available.  However, 
attempts  we  have  made  to  find  a common  correlation  of  these  data  for 
example  by  a single  turbulence  model,  have  not  been  successful.  The 
problem  is  complicated  by  very  complex  shock  wave  interactions  in  the 
near  flowfield.  However,  even  a generalized  description  of  the  far  flow- 
field  seem.s  unavailable  at  this  time. 

Available  experimental  data  for  two  phase  jets  includes  very  limited 
studies  on  incompressible  (9)  and  compressible  (ID)  jets.  These  data 
are  generally  not  sufficient  to  verify  analytical  descriptions  of  the 
flow  since  only  a limited  range  of  conditions  have  been  studied.  None- 
theless, analytical  two  phase  flow  solutions  have  been  obtained  by  limited 


solation  of  the  Navier-Stokes  Equations  (3,  5,  11).  In  order  to  describe 
these  flows  more  accurately,  experimental  measurements  in  two  phase  axi- 
symmetric  jets  must  be  made  over  a wider  range  of  experimental  conditions. 

The  effects  of  particles  in  real  gas  flows  are  reflected  in  a 
number  of  changes  in  the  flow.  In  general,  particle  velocity  and  tern-  j 

peratJi'e  may  be  different  from  the  gas  velocity  and  temperature,  the 
particles  may  diffuse  at  different  rates  in  shear  layers,  and,  finally, 
m.ass  may  be  exchanged  between  the  particle  and  the  gas.  Attempts  to 
account  for  these  effects  have  led  to  a number  of  approximate  descrip- 
tions of  the  flow.  Particle  velocity  and  temperature  have  been  assumed 
to  be  in  equilibrium  with  the  gas  or  frozen  at  its  initial  velocity  or 
temperature.  The  same  equilibrium  or  frozen  description  has  been  applied 
to  particle  diffusion  effects  and  to  phase  change  or  chemical  reaction  of 
the  particles.  These  two  limiting  conditions  are  most  useful  as  extremes 
of  the  real,  dynam.ic  interactions  between  the  two  phases. 

In  studies  under  this  grant,  an  analytical  solution  procedure  was 
developed  for  the  description  of  particle  flows  in  shear  layers.  The 
model  incorporates  the  dynamics  of  interaction  between  particle  and  gas 
phase.  Governing  conservation  equations  for  gas  and  solid  phases  are 
solved  by  a modification  of  the  integration  procedure  of  Patankar  and 
Spalding  (ll"!  for  2D  parabolic  flow.  Analytical  results  were  compared 
to  detailed  flow  and  particle  concentration  measurements  made  for  jet 
flow  and  for  duct  flow.  It  was  found  that  the  analytical  mean  solid 
fluxes  in  a two-phase  jet  flow  agree  well  with  the  experimental  result'; 
and  that  the  analytical  velocity  profiles  in  a two-phase  duct  flow  agr(;c 
fairly  well  with  the  experimental  results.  The  stiidy  also  considered 
overall  properties  of  the  jet  flowfields  including  spreading  rates,  en- 
trainment formulae  and  centerline  velocity  decay. 

hi  
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4. 

final  investigation  reported  here  is  a study  of  the  chemistry 
of  afterburning  plumes.  Chemical  reactions  in  axisymmetric  jets  similar 
to  jet  and  rocket  exhaust  arise  from  mixing  of  unburned  fuel  in  the  jet 
vvith  air.  The  complete  description  of  the  flowfield  requires  knowledge 
• of  viscosity  and  chemical  kinetics  in  the  jet.  Recent  measurements  of 

plume  radiation  properties  (13)  and  predictions  based  on  calculations  of 
reacting  flowfields  have  not  been  in  agreement.  One  possible  explanation 
of  this  discrepancy  has  been  the  possibility  of  a mixture  at  the  nozzle 
exit  with  a substantial  e/.cess  of  unburned  fuel.  Calculations  which  we 
made  for  reacting  jets  with  a nonequilibrium  nozzle  exit  composition 
qualitatively  supported  this  explanation  for  the  radiation  discrepancy 
but  could  not  quantitatively  account  for  the  difference. 
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II.  Turbalence  Modeling  of  Shear  Layers 


The  problem  of  turbulent  mixing  of  compressible  jets  with  air  is 
complex  and  many  practical  data  are  not  available.  Although  consider- 
able study  has  been  undertaken  by  other  authors  of  incompressible  jets  (25) 
both  experimental  and  analytical  study  of  compressible  jets  is  much  less 
complete.  An  historical  review  of  turbulence  models  for  axisymmetric 
jets  is  presented  in  Table  I.  Some  of  these  results  will  be  discussed 
here  with  relevance  to  compressible  afterburning  jet  plumes. 

Harsha  (23)  has  critically  reviewed  many  of  the  turbulence  models 
studied  so  far  an^  recommends  that  attempts  to  modify  basic  Prerdtl 
eddv  viscosity  model  or  the  mixing  length  theory  to  make  them  apply  to 
more  complex  flows  is  not  productive.  None  of  the  modifications  of  the 
Prer.dtl  eddy  viscosity  model,  including  Donaldson  and  Grey  (7)  compress- 
ibility correction  are  capable  of  greatly  altering  the  basic  shape  of  the 
axial  centerline  velocity  decay  curve,  and  the  shape  predicted  by  the 
Prendtl  model  and  all  of  its  derivatives  is  incorrect  for  complex 
(twc-gas)  flows.  On  the  ether  hand,  the  displacement-thickness  model 
proposed  by  Schetz  (21)  is  the  only  locally-dependent  model  to  show  the 
proper  behavioral  trends  for  hydrogen-air  mixing.  Because  of  this,  its 
use  should  be  investigated  in  other  dissimilar-gas  flows.  Given  some 
knowledge  of  the  initial  turbulent  shear  stress,  the  turbulent  kinetic 
energy  method  (24)  is  capable  of  providing  better  and  more  uniform  pre- 
dictions over  a wider  range  of  flows  than  any  other  models  investigated. 
Because  of  this,  it  clearly  holds  the  greatest  promise  for  future  devel- 
opment . 

In  order  to  evaluate  the  predictive  ability  of  some  of  these  tur- 
bulence models,  calculations  have  been  made  for  conditions  representative 


Author/McxJel  Year  i.'xpression  for  Eddy  Viscosify  Remarks 
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cl  some  available  data  for  compressible  jets.  These  calculations  were 
made  with  the  solution  procedure  of  Reference  14  with  several  turbulence 
.models.  Previously  the  program  was  written  with  an  explicit  linite 
difference  scheme  but  it  was  later  modified  to  use  an  implicit/explicit 
finite  difference  scheme  which  resulted  in  more  efficient  computation. 

The  progra.m  predicts  gas  dynamic,  chemical  and  electrical  properties  of 
axisymjr.etric  rocket  plumes  from  sea  level  through  the  continuum  flow 
regime . 

Centerline  velocity  decay,  spreading  rate  of  jet  and  species  con- 
.er.tratior  have  been  calculated.  Some  experimental  results  are  avail- 
able in  references  (23)  and  (25,  26,6,  7,8)  for  compressible  jet-into- 
still  Air,  coaxial  Air-Air  mixing,  coaxial  Hydrogen-Air  mixing  (11), 
heated  air  and  cooled  air  (6),  Nitrogen,  Methane  (7)  with  which  we  f-an 
compare  our  results  (Figure  1).  While  the  nature  of  variation  of  center- 
line  velocity  looks  in  order,  it  does  not  give  the  expected  type  of 
variation.  From,  the  experimental  results  available  it  is  anticipated 
that  centerline  velocity  should  decay  by  about  70%  at  around  40  ratHi 
but  it  is  not  happening  in  our  case  with  any  of  the  viscosity  models  used. 
The  spreading  rate  is  also  not  quite  of  the  form  as  given  by  Warren  (f>). 
This  further  suggests  the  need  of  experimental  as  well  as  theoretical 
study. 
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III.  Particle  Flows  in  Shear  hayerr, 

The  presence  of  particles  in  exhaust  plumes,  particularly  turbulent 
mixing  regions,  can  contribute  to  design  and  prediction  problems.  Par- 
ticles in  the  exhaust  can  contribute  band  radiation  in  excess  of  that 
characteristic  of  molecular  radiation.  The  region  of  plume  separation 
from  the  vehicle  can  further  act  as  a flame-holder  for  the  afterburning 
plume.  This  flameholding  is  sustained  by  shear  layer  reactions,  including  | 

oxidation  of  carbon  particles,  rather  than  the  recirculation  zc!ie  created  ] 

i 

by  the  plume  separation.  The  accurate  overall  description  of  the  flow  j 

and  reaction  phenomena  must  consider  particle/gas  flow  and  reaction.  ; 

Major  unresolved  difficulties  persist  in  such  a description.  In  addition 
to  particle/gas  mixing  phenomena  which  determine  initial  temperature  and 
velocity  profiles  at  the  nozzle  exit,  the  description  of  particle  forma- 
tion and  particle/gas  reactions  are  little  understood.  The  work  wftic:h 
is  reported  here  is  an  approach  at  analytical  description  of  some  of  < 

these  phenomena.  | 

The  effects  of  particles  in  real  gas  flows  are  reflected  in  a number  | 

of  changes  in  the  flow.  In  general,  particle  velocity  and  temperature  may  | 

be  different  from  the  gas  velocity  and  temperature,  the  particles  may 

I 

diffuse  at  different  rates  in  shear  layers,  and,  finally,  mass  may  be  | 

exchanged  between  the  particle  and  the  gas.  Attempts  to  account  for 
these  effects  have  led  to  a number  of  approximate  descriptions  of  the 
flow,  Table  II.  Particle  velocity  and  temperature  have  been  assumed  to  ! 

I 

be  in  equililrium  with  the  gas  or  frozen  at  its  initial  velocity  or  tem- 
perature. The  same  equilibrium  or  frozen  description  has  been  af)|-lie(i 
to  particle  diffusion  effects  and  to  phase  change  or  chemical  reat  l ion 
of  the  particles.  These  two  limiting  conditions  are  most  useful  as  ex- 
tremes of  tlie  real,  dynamic  interactions  between  the  two  piiases. 


Table  II,  Particle  properties  relative  to  gas  in 
two  phase  nozzle  and  exhaust  flow. 


EQUILIBRIUM 

KINETIC 

FROZEN 

VELOCITY 

B 

A,C,D,E 

B 

TEMPER. 

B 

C,D,E 

B 

DIFFUSION 

G 

A 

PHASE  CHG. 

G 

E.G 

B Altman  and  Carter,  1956 
C Kliegel  and  Nickerson,  1967 
D Crowe  and  Pratt,  1972 
E Jensen  and  Wilson,  1974 

F Pergament,  1 974 

G Genovese  et  al . , 1971 


Table  m.  Characteristic  data  for  exhaust  plumes. 


Gas  Temperature,  Tg 

800  - 2500“ K 

Gas  Pressure,  P 

0.1  - 1.0  atm 

Gas  Velocity,  Ug 

2.5  km/sec 

Gas  Density,  pg 

2 - 6 X lO”**  gm/cc 

Particle  Denisty,  Pp 

2.0  gm/cc 

Particle  Size,  dp 

.1  - 10  ym 

Particle  Loading 

5 - 25% 

Particle  Re 

10“^  - 10-**  lu  - u 1 

9 P 

Temperature  Gradient 

10®  - 10®  “K/sec 

Velocity  Gradient 

10®  - 10®  cm/sec/sec 

Transit  Time 

10-“  - 10'®  sec 

! 
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A . Conservation  Equations  for  a Multiphase  Flow 

The  conservation  equations  for  two-dimensional  particle  laden 
boundary  layers  can  be  written  as  follows  (33): 

Conservation  of  mass  for  carrier  fluid  including  chemical  reaction: 


c^PfUf 


aPfVf 

~§y~ 


= r 


Cl) 


where  r is  the  rate  of  conversion  of  particles  to  gas, 
Conservation  of  mass  for  solid  particles: 

3p  u dp  V 

-£-2  + P P =_r 

?x  ^ dy  ‘ 


(2) 


Conservation  of  momentum  for  carrier  fluid  (neglecting  the  momentum 
due  to  generated  gas): 

d 


du^  du. 

f f 

^x  f dy 


,,2 

d Uf  P 

U f p 
+ T + n + 

dx  ^f 


(3) 


m 


where  = particle  velocity  relaxation  time  = Ppd^"^/18|j.  . 


du  du  u -Ur 

u _E  + V — £ = £ ^ 

p dx  p dy 


Conservation  of  momentum  in  the  axial  direction  for  the  particles: 

(4) 

Conservation  of  momentum  in  tlie  normal  direction  for  the  particles; 

(5) 


dV  'V  V -Vr 

U ^ + v^  = -P.  1 

P ax  p jy  X 

^ m 


Since  the  momentum  transferred  in  the  normal  direction  is  small 
compared  to  that  in  the  axial  direc.ion,  Equation  (5)  will  be  neguected 
in  the  following  analysis. 


To  describe  analytically  a suspension  flow  in  a viscous  layer, 
filiations  (3)  and  (d"),  wtiich  are  coupled,  a long  with  Equations  ( 1;  and  (2)  , 
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which  are  also  coupled,  have  to  be  solved  simultaneously.  A relation- 
ship exists  between  u^  and  u^  , assuming  that  particle  inteiactions 
are  negligible  and  the  relative  turbulent  drag  is  small,  e.g.  (33): 

.15g-^V*^\p  -pp*^ 

(u^-u  ) = [ 3 (€ 

Pf 

Where:  d^  = diameter  of  solid  particles 

Dp  = density  of  solid  phase 
0^  = density  of  gas  phase 

Substitution  of  this  equation  into  Equation  (3)  yields: 


""f  "57  ^ "'f  “§7  - (^3 

This  equation  differs  from  the  equation  of  conservation  of  axial  momenti 
for  single  phase  flow  only  in  the  source  term  d where: 

d = - -i 

pu  dx 

for  single  phase  flow,  and 


u^  L p^  dx  p^ 


,15g*^^d^‘^^(p  -p.)'^ 
^ p ^ ^p  f 


,.29,, .43 


“f 


3^} 


for  two- phase  flow. 

This  observation  suggests  that  one  way  to  solve  the  two- phase  flow 
problem  would  be  by  modifying  the  source  term  in  the  solution  procedures 
for  one- phase  flow  and  solving  for  the  velocity  of  gas,  then  using  equa- 
tion (1)  and  solving  for  the  velocity  of  the  particle  phase.  The  solution 
can  be  obtained  by  a modification  of  the  integration  algorithm  developed 
by  Patankar  and  Spalding  (12)  for  2D  parabolic  flow. 
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B . Calculated  Results 

The  two-phase  jet  problem  has  been  treated  experimentally  by 
Lewis  (34)  who  measured  axial  and  radial  particle  flow  profiles  at 
particle  Reynolds  numbers  in  the  transition  regime.  Data  were  correlatoti 
in  terms  of  decay  of  centerline  particle  flux,  i.e., 


(uc) 

(uc)' 


c-t 


o 

n 

o TTC  X 

and  radial  profiles  of  particle  flux 

(uc). 


(UC)^  . 


cl 

The  mean  solid  fluxes  Uc  above  are  expressed  in  terms  of  an  empirical 
constant  c = .0596  and  cross-sectional  area  A at  different  axial  (x) 
and  radial  (r)  locations. 

As  an  illustration  of  the  present  analysis,  calculations  were  made 
for  the  gas- solid  mixture  in  Lewis'  experiments  and  the  required  para- 
meters needed  for  the  Patankar- Spalding  program  were  evaluated  on  that 
basis  for  different  particle  loading  ratios.  The  program  was  modified 
for  a jet  and  the  initial  conditions  and  the  calculated  properties  of  the 
gas/particle  mixture  were  incorporated  into  it.  The  value  used  for  the 
viscosity  of  the  gas  was  that  of  air  and  far  from  the  center- line  of  the 
jet  it  was  set  equal  to  twice  that  of  air,  because  of  higher  rate  of 
momentum  transfer  between  the  two  peases  in  that  region.  Figures  2 and 
3 show  the  results  for  different  loading  ratios. 

The  two-phase  duct  flow  problem  was  treated  experimentally  by  Doig 
and  Roper  (35),  who  measured  the  velocity  of  air  as  a function  of  radial 
viistance  for  tiifferent  axial  distances  in  a duct  of  1.7  inch  insirle  diam- 
eter. Splierical  particles  were  used  with  loading  ratio  ranging  from  .1 
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Figure  2.  Calculated  and  experimental  particle 
flux  in  a spreading  two  phase  jet.  tx 
perimental  data  by  Lewis  134;  for  low 
particle  loading. 


Figure  3.  Calculated  and  experimental  particle 
flux  in  a spreading  two  Ph^se Jet. 
Experimental  data  by  LewisV34J  for 
moderate  particle  loading. 
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to  5.5.  Data  were  plotted  in  terms  of  velocity  versus  square  of  rela- 

2 

tive  distance  from  axis  (2r/D)  . The  properties  of  the  flow  in  the 
experiments  were  calculated  for  different  loading  ratios,  Figures  4 and 
5,  for  the  turbulent  transport  properties  determined  by  the  jet  flow 
calculation.  In  these  calculations  a mixing  length  turbulent  viscosity 
was  used  for  the  gas: 


u 


t 


5u 

?y 


where  I = 0.09  6 and  5 is  the  jet  half  width. 

This  rather  simplified  model  for  the  two-phase  jet  and  duct  flow 
is  seen  to  give  generally  close  agreement  to  the  experimental  data.  The 
results  indicate  the  feasibility  of  the  analysis  as  an  approximation  to 
two- phase  flow  in  a parabolic  boundary  layer  region  of  an  exhaust  flow- 


field. 


r 


FROM  AX1S,(^)^ 

Figure  4.  Velocity  profile  in  core  region  of  two 
phase  duct  flow  with  loading  ratio  of  3. 
Ref.  (35). 


FROM  AXIS, 


Figure  5.  Velocity  profile  in  core  region  of  two 
phase  duct  flow  with  loading  ratio  of 
4.4.  Ref.  (35). 


; 


?0. 

IV.  Chemical  No.nequilibrium  Effects  in  Shear  Layers 

Recent  measurements  by  Ebeoglu,  et  £l.  (13)  have  indicated  sig- 
nificant discrepancies  between  measured  plume  radiation  characteristics 
and  predictions  based  on  the  methods  outlined  above.  Measurements  of 
radiant  intensity,  spatial  distribution  of  radiant  intensity  and  plume 
radiation  spectral  distribution  were  made  on  a small  kerosene/cxygen 
rocket.  The  measurements  were  then  compared  to  predictions  using  a 
Naval  Weapons  Center  plume  prediction  code  and  an  Air  Force  Armament 
Laboratory  code.  It  was  found  that  neither  code  predicted  the  plume 
radiation  characteristics  with  reasonable  accuracy.  At  oxidizer/fuel 
ratios  less  than  three  the  measured  radiant  intensity  in  the  4 tc  4.5  (j,m 
band  was  as  much  as  an  order  of  magnitude  greater  than  predicted  (Figure  6). 
For  0/F  ratios  greater  than  4 the  calculated  radiation  was  twice  that 
measured.  Measurements  of  spectral  distribution  showed  bands  centered 
at  2.7  and  4,3  pm  in  agreement  with  predictions  which  considered  CO^  and 
H^O  as  the  dominant  radiating  species.  However,  the  calculations  predict 
less  radiant  intensity  in  the  4.3  pm  band  as  compared  to  the  2.7  pm  region. 
The  measurement  of  spatial  distribution  of  radiant  intensity  also  showed 
disagreement  between  calculations  and  experiment.  The  peak  radiant  in- 
tensity was  predicted  to  be  much  closer  to  the  nozzle  exit  than  was 
measured  from  the  model  rocket  motor. 

These  resu.lrs  suggest  serious  deficiencies  in  the  ability  of  state 
of  the  art  computer  codes  to  predict  plume  radiation  properties.  The 
discrepancies  may  arise  from  several  calculational  or  experimental  dif- 
ficulties. The  analysis  above  assumed  combustion  to  chemical  equilibrium 
products  prior  to  expansion  in  the  nozzle.  Since  the  plume  properties 
at  the  exit  plane  were  not  measured,  the  initial  conditions  for  the 
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Figure  6.  Variation  of  Radiant  Intensity  with  0/F  Ratio  (Ref.  13) 


nr 


afterburning  plume  were  not  known  and  should  be  determined.  Separate 
calculations  of  nonequilibrium  nozzle  flow  by  Vamos , £t  a^.  (39)  sug- 
gested that  if  equilibrium  combustion  is  assumed,  the  nozzle  exit  con- 
centration of  CO^  and  HgO  are  below  equilibrium.  Predictions  based  on 
these  results  would  indicate  lower  initial  radiation  intensity.  The 
possibility  of  an  increased  contribution  to  radiation  intensity  due  to 
afterburning  in  excess  of  that  calculated  has  been  suggested  (13). 

In  an  attempt  to  describe  the  plume  chemical  characteristics  and 
hence  more  accurately  predict  radiation  properties,  this  study  has 
focused  on  the  exhaust  plume  initial  composition.  Calculation  of  the 
composition  history  in  the  exhaust  plume  has  been  made  using  the  avail- 
able chemical  kinetic  data  and  considering  two  composition  initial 
profiles.  In  the  first  case  the  composition  begins  at  the  nozzle  exit 
equilibrium  composition  and  reacts  as  air  is  entrained  into  the  jet. 

In  the  second  case  the  composition  at  the  nozzle  exit  is  adjusted  to 
simulate  an  excess  of  unburned  fuel. 


A . Calculated  Results 

The  NASA  CEC  Computer  program  (36)  was  used  to  calculate  the 
equilibrium  combustion  products  in  the  chamber  and  exit  plane.  Calcula- 
tions were  made  for  a chamber  pressure  of  100  psia  and  initial  temperature 
of  300 °k.  The  principal  constituents  and  their  mole  fractions  are  shown 
in  Figure  7.  Other  minor  species  not  shown  in  this  figure,  such  as  HCO, 
HO2,  0 and  H were  predicted  and  included  in  the  kinetic  calculations. 

HCO  is  important  in  rich  0/F  ratios,  while  HO2  appears  only  in  the  lean 
0/F  ratio  mixtures. 


F RATIO 


1 composition  in  the 
vs  . chamber  0/F  ratio. 
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The  fast  computer  program  for  nonequilibrium  rocket  plume  predic- 
tions CD  has  beer,  used  to  compute  the  thermodynamic  and  chemical 
properties  of  a rocket  exhaust  plume.  The  rocket  performance  input  data 
and  chemical  species  concentration  were  given  by  the  equilibrium  calcula- 
tions. The  chemical  kinetic  data  were  obtained  from  references  (2),  (4) 
and  (5)  and  are  reproduced  in  Table  IV.  The  same  engine  characteristics 
as  those  of  (2)  were  used  in  these  calculations.  ; 

The  nonequilibrium  chemical  species  concentration  data  used  in  : 

this  analysis  were  obtained  by  manipulating  the  equilibrium  data.  For  1 

example,  20%  nonequilibrium  condition  was  achieved  by  20%  increase  in  j 

CO,  O2  and  mole  fraction  ancf  20%  decrease  of  CO2  and  H2O  from  the 
equilibrium  data  obtained  by  NASA  computer  program.  [Figure  7]  The  | 

validity  or  accuracy  of  this  assumption  is  limited,  but  qualitatively  | 

it  follows  the  existing  data.  | 

1 

I 

I 

B . Discussion  | 

1 

The  primary  results  of  the  calculation  are  presented  in  Figures  8 \ 

I 

and  9.  In  Figure  8 the  axial  variation  of  centerline  concentration  of  ^ 

CO2  is  plotted  for  several  overall  0/F  ratios  beginning  with  an  equilib-  \ 

rium  composition  at  the  nozzle  exit.  At  the  most  fuel  rich  mixture  condi- 
tion CO2  increases  very  rapidly  with  axial  distance  which  is  similar  to 
measurement  of  total  radiation  by  Ebeoglu,  et  al.  Concentration  at  low 
0/F  was  also  found  to  exceed  values  at  higher  0/F  by  as  much  as  a factor 
of  two.  If  infrared  radiation  is  taken  to  be  primarily  due  to  CO2  and  i 

K2O,  this  suggests  an  increase  in  radiation  at  low  0/F  ratios.  Calcula- 
tions for  radial  profiles  comparing  the  equilibrium  and  nonequilibrium 
initial  condition  are  given  in  Figure  9,  These  data  show  a further  in- 
crease in  CO2  ho  the  nonequilibrium  initial  condition. 


REACTIONS  BEING  CONSIDERED  KR:=A*EXP(  B/ RT)/T 
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Figure  8.  Axial  variation  of  CO^  concentration  in 
axisymnietric  jet  plume  at  various  0/F  ratios. 
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